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1.  INTRODUCTION 


J 

The  near-millimeter  region  of  the  spectrum  (NMMW:  0.1-1  THz)  is  an  active 
area  for  research.  Possible  applications  lie  in  short-range  communications, 
radar,  radiometry,  and  radio  astronomy.  Atmospheric  effects  of  transmission 
and  emission  dominate  performance  simulations  and  the  complex  refractivity  N 
provides  a  measure  of  the  interactions  between  radiation  and  the  atmospheric 
propagation  medium.  A  reliable  N-model  allows  calculation  of  frequency-depen¬ 
dent  rates  for  delay  (real  part)  and  attenuation  (imaginary  part)  based  on 
measurable  meteorological  variables.  Dry  air  and  atmospheric  water  vapor  are 
major  millimeter-wave  absorbers;  so  are  suspended  droplets  (haze,  fog,  cloud) 
and  precipitating  water  drops  that  emanate  from  the  vapor  phase.  Laboratory 
research  and  analytical  studies  have  been  conducted  with  the  primary  purpose 
of  understanding  NMMW  attenuation  |cu-db/km) c  and  delay (B^_ps/km)c  rates. ^ 
Emphasis  was  placed  on  the  fundamental  concepts  that  support 
an  a/0  formulation. 

Refractivity  N  for  moist  air  can  be  obtained,  in  principle,  by  a  line-by¬ 
line  summation  over  all  molecular  absorption  lines.  In  practice,  various 
degrees  of  approximations  are  employed  to  reduce  labor  and  computer  time 
required,  since  the  number  of  contributing  spectral  lines  by  the  dominant 
absorbers  (water  vapor  and  oxygen)  and  by  various  trace  gases  (e.g.,  O3) 
exceeds  10,000.  A  practical  NMMW  model,  indexed  MPM,  consists  of  local  H20 
(30x)  and  02  (48x)  lines  below  1  THz  and  an  approximation  to  the  contributions 
by  H20  lines  above  1  THz  [3],  [6],  [8],  [9],  [Appendix]. 

The  experiments  have  been  performed  at  138  GHz  to  measure  absolute  atten¬ 
uation  rates  by  dry  air,  moist  air,  water  vapor,  and  water-vapor  mixtures  with 
nitrogen  (N2),  oxygen  (02),  and  argon  (Ar)  at  temperatures  between  8  and  43°C, 
total  pressures  up  to  1  atm  and  higher,  and  relative  humidities  up  to  95  per¬ 
cent.  The  experimentally  observed  absorption  is  not  described  by  standard 
line  shape  models.  Such  failure  reveals  difficulties  in  modeling  frequency, 
temperature  and  pressure  dependencies  for  moist  air  attenuation.  An  unex¬ 
plained  excess  is  identified  for  which  the  name  "water-vapor  continuum"  was 
coined  since  it  appears  to  increase  smoothly  with  frequency  within  the  NMMW 


Experimental  studies  inspired  model  formulations  along  practical  lines. 
The  MPM  program  is  a  user-fr iendly,  PC -operated  code  that  generates  numerical 
values  of  a(f)  and  8(f)  for  frequencies  f  up  to  1000  GHz.  Input  parameters 
are  five  measurable  atmospheric  quantities:  barometric  pressure  P,  ambient 
temperature  T,  relative  humidity  RH,  suspended  droplet  water  content  w,  and 
rainfall  rate  R.  Controlled  laboratory  measurements  were  limited  to  moist  air 
studies  (P,  T,  RH) ,  and  the  data  obtained  at  138  GHz  are  reasonably  complete 
and  accurate  to  assess  water  vapor  pressure  and  temperature  dependencies  for 
the  water-vapor  continuum.  Both  variabilities  point  to  the  distinct  poss¬ 
ibility  of  an  absorption  mechanism  related  to  water  vapor  but  not  accounted 
for  by  molecular  theory  of  H2O. 

This  report  is  organized  in  three  parts.  The  first  part  (Section  2) 

gives  details  of  the  experimental  setup,  its  achieved  performance,  and  a 

summary  of  reduced  data.  After  many  improvements,  a  detection  sensitivity  of 

a  .  «  0.05  dB/km  or  1.2  x  10“7cm~*  was  realized.  In  the  second  part  (Sec- 
min 

tion  3),  results  from  the  laboratory  experiments  are  applied  (a)  to  calibrate 
the  MPM  program  with  an  empirical  continuum  term,(b)  to  demonstrate  the  para¬ 
metric  flexibility  of  the  code  ( i . e . ;  f,  RH,  and  P  can  be  selected  as  vari¬ 
ables),  and  (c)  to  conjecture  on  the  physical  basis  of  a  NMMW  water  vapor 
continuum  that  is  defined  by  the  limited  H20  line  base  of  MPM.  A  final  sec¬ 
tion  4  contains  examples  of  recently  reported  data  from  laboratory  and  field 
experiments  on  water  vapor  absorption  (10-430  GHz)  and  their  comparison  with 
MPM  predictions.  Overall  agreement  is  good,  particularly  in  atmospheric 
transmission  windows. 
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2.  LABORATORY  STUDIES  OF  MOIST  AIR  ABSORPTION  AT  138  GHZ 


Controlled  experiments  that  simulate  atmospheric  conditions  provide  test 
cases  for  studying  specific  contributions  to  N  in  isolation.  Assessments  of 
basic  physical  principles  underlying  the  attenuation  rate  o  are  difficult  to 
make  from  measurements  in  the  actual  atmosphere.  The  objective  of  this  study 
was  to  measure  water  vapor  continuum  absorption.  A  test  frequency  of  138  GHz 
was  selected  because  of  its  remoteness  from  local  H20  lines.  The  expected 
window  attenuation  falls  in  the  range  0.1  to  5  dB/km  and  the  required  detec¬ 
tion  sensitivity  calls  for  a  long  ( >0 . 1  km)  effective  path  length,  which  can 
be  attained  with  a  resonant  absorption  cell. 

The  response  curve  A(f)  of  an  isolated,  high  Q-value  resonance  is  detec¬ 
ted  with  a  power  (square-law)  detector.  Both,  the  peak  value  a0  at  center 
frequency  f p  and  the  bandwidth  b0  spread  over  a  range  fp  ±  bQ/2  at  the  level 
a0/2  might  be  used  to  detect  the  relative  attenuation, 

ar  -  8.686(/a^7a  -1)  -  8.686 (b/bQ  -1)  dB,  (1) 

of  an  absorbing  gas,  that  changes  the  corresponding  quantities  to  a  and  b  when 
introduced  into  the  resonator.  Around  138  GHz  it  is  possible  to  design  a 
compact  (20  cm  mirror  spacing)  Fabry-Perot  resonator  with  a  loaded  Q-value  on 
the  order  of  4  x  105,  that  defines  (Q  -  fR/bQ)  a  resonance  bandwidth  of 
bo^350  kHz. 

A  crucial  question  to  be  resolved  is  if  amplitude  (aQ/a)  or  frequency 
(b/b0)  detection  schemes  provide  the  optimum  sensitivity  for  the 
spectrometer.  After  exhaustive  testing  it  was  found  that  digital  averaging  of 
A ( f )  ,  displayed  over  a  frequency  span  AfM  ■  fR  t  6bQ  was  capable  of  resolving 
a0/a^  1.002  (512  pts)  -  but  only  b/bQ^  1.015  with  1024  pts  .  Amplitude 
peak-value  detection  provided  optimum  sensitivity  for  absorption  studies.  The 
frequency  span  Af^  is  needed  to  establish  the  baseline  A(f)  -  0  of  the  reso¬ 
nance  response.  In  addition,  a  detection  at  fp  eliminates  corrections  for 
dispersive  distortions  of  A(f)  [21], 

Absolute  calibration  of  absorption  is  accomplished  by  defining  an  equi¬ 
valent  path  length  (b0  in  kHz), 

Le  -  47.71/b0  km,  (2) 
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for  the  resonance  spectrometer  operating  at  fR.  From  (1)  and  (2)  follows  that 
the  absolute  power  attenuation  rate  of  an  absorbing  gas  is  given  by 

a  -  0.1820bo(/a^7a  -  1)  dB/km.  (3) 

If  the  assumptions  made,  so  far,  can  be  brought  to  fruition,  then  the  pro¬ 
jected  detection  sensitivity  (aQ/a  -  1.002,  b0  -  350  kHz)  is 
“min  *  0*064  dB/km,  which  is  adequate  for  che  planned  water  vapor  studies. 

2.1  Experimental  Arrangement 

The  measuring  system  consists  of  the  millimeter  wave  resonance  spectro¬ 
meter  and  a  humidity  simulator.  An  insulated  box  contains  a  high-vacuum 
stainless  steel  vessel  which  houses  a  temperature-controlled  mini-lake  (10  cm 
across)  and  the  resonator.  Temperatures  are  controlled  to  better  than  1/100 
of  a  degree  Celsius;  pressure  ranges  over  seven  orders  of  magnitude  (102  to 
10-5  kPa);  and  relative  humidity  is  varied  between  0  and  99.5  percent. 

Schematical  diagrams  of  physical  and  electronic  (Figure  1)  arrangements 
convey  an  overview  of  the  experiment.  A  temperature-controlled  water  reser¬ 
voir  serves  as  the  vapor  source.  Electropolished  stainless  steel  was  used 
exclusively  as  construction  material.  Various  hydrophobic  coatings  were 
studied  as  possible  means  for  neutralizing  the  absorption/desorption  cycle  of 
surfaces  exposed  to  water  vapor  [1],  but  were  abandoned  in  favor  of  slightly 
heating  the  mirrors  of  the  resonator.  Four  fast-responding  (t  <  Is)  temper¬ 
ature  sensors  inside  the  cell  diagnose  any  disturbance  of  the  gas  condition. 
Data  acquisition  was  computer  controlled. 

The  resonator  inside  a  vacuum  chamber  (Figure  lc)  is  the  heart  of  the 
absorption  spectrometer .  A  key-word  summary  of  its  specifications  reads  as 
follows:  Fabry-Perot  reflection-type,  semi-confocal  arrangement,  10  cm  mirror 
diameter,  mirror  heating:  l°C/0.3  W;  Fresnel  number:  6;  coupling  factor 
k  -  0.0550,  coupling  hole:  circular,  0.65  mu  diameter,  0.075  mm  double  Mylar 
vacuum/pressure  seal;  resonance  frequency,  selected  for  optimum  performance  of 
the  available  klystron,  fR  -  137.80  GHz,  stability  Af /f R-  1.1  x  10-6/'C 
temperature  compensated  and  insensitive  to  pressure  loads  from  0  to  2  atm; 
resonance  bandwidth:  b0  -  334  kHz  yielding  an  effective  path  length, 
t_£  -  0.141  km  (Eq.  2);  mirror  spacing:  x^  -  182  A/2  -  198.0  mm;  micrometer 


tuning:  0.3175  mm/rev  ■  2500  SU/rev  with  a  resolution,  AxR  -  0.127  pm/SU, 
which  converts  into  a  frequency  change, 

AfR  -  -fR(Ax/x)R  -  -88.4  kHz/SU.  (4) 

Over  a  2.5  cm  tuning  range  of  the  micrometer,  the  resonance  at  xR  was  the 
"best”  out  of  52  choices  with  respect  to  other-mode  interferences  and  Q- 
values. 

The  detector  of  the  resonance  pulse  A(t)  was  an  InSb  bolometer,  cooled  to 
4.2  K-LHe,  and  connected  by  a  WR-12  waveguide  mount.  The  maximum  voltage 
output  of  the  preamplifier  was  10  V.  A  power-linear  response  of  4.2  V/mW  was 
measured  up  to  1.0  V.  With  a  50  kHz  detection  bandwidth, the  noise  power  was 
about  5  nW.  The  bolometer  bias  (92  mV)  served  as  cryogenic  thermometer 

(T0). 

The  power  source  was  a  138  GHz  klystron  (20  mW)  that  was  frequency- 
modulated  by  a  sawtooth  voltage  generator  to  provide  frequency-to-time  domain 
conversion.  The  modulation  frequency  was  exactly  500  Hz  derived  from  a  10  MHz 
frequency  standard.  The  linear  sawtooth  ramp  was  gated  at  exactly  1925  ps, 
the  fly-back  time  took  only  2  ps,  which  eliminated  A(t)  from  the  retrace.  The 
modulation  sensitivity  of  the  klystron  was  determined  accurately  to  be 
12.56  MHz/V  by  using  the  resonance  peak  a0  as  frequency  marker,  mechanically 
tuned  to  the  end  points  of  the  ramp  and  checked  for  linearity  over  a 
modulation  voltage  range  from  0.500  to  3.000  V.  The  tuning  uncertainty  of 
0.5  SU  introduced  about  0.5  percent  error  in  the  frequency-to-voltage 
conversion  factor  (12.56  MHz/V). 

The  reference  power  aQ  -  k-aK  for  (1)  was  measured  by  periodically  (30s) 
switching  on  computer  command  to  a  36  V  modulation  voltage,  which  displayed 
the  power  mode  AK(t)  with  a  peak  value,  aK  -  909  mV  and  a  bandwidth, 
bK  -  186  MHz.  This  feature  provided  an  automatic  calibration  of  ap  (1)  when 
the  power  level  changed  during  a  data  run  (e.g.,  refractive-tun ing  reduces  fR 
and  the  peak  value  a^  can  change  when  the  klystron  is  readjusted). 

An  electronic  lock-on  circuit  kept  the  rsonance  fR  centered  with- 
time  frame;  that  is,  the  klystron  center  frequency  f^a^)  was  prevented  from 
drifting  with  reference  to  fR-  A  baseline  A(t)  -  0  was  established  by 
adjusting  fK  in  such  a  way  that  two  adjacent  A(t)  frames  were  displayed  on  a 
control  scope  and  balanced.  The  reflected  resonance  signal  A ( t )  could  be 


eliminated  by  injecting  a  steel  rod  into  the  resonance  volume.  This  measure 
provided  baseline  reference  data  in  the  "low"  (0.36  V)  and  more  accurate  peak 
readings  aK  in  the  "high"  (36  V)  mode  of  the  modulator. 

The  waveform  processor  for  A(t)  was  a  digital  storage  oscilloscope, 
synchronized  with  the  sawtooth  modulator  and  capable  of  resolving  512/1024  pts 
per  2000  ps.  Operational  resolution  was  typically  3.91  ps/pt.  The  modulator 
voltage  for  A(t)-detection  was  0.360  V  resulting  in  a  frequency  resolution  of 
9.2  kHz/pt,  which  is  an  improvement  over  (4).  Extensive  averaging  of  the 
repetitive  waveforms  A ( t )  and  A^(t)  up  to  200  times  was  performed  in  real  time 
to  improve  the  S/N  of  the  a  and  b  results. 

For  a  measurement  of  the  resonance  bandwidth  bQ,  the  digitizing  incre¬ 
ments  were  doubled  (1024  pts)  and  the  modulation  voltage  was  varied  in  the 
0.150  to  0.400  V  range,  allowing  the  resolution  uncertainty  to  be  reduced  to 
±5  kHz.  The  resulting  error  in  an  absolute  calibration  of  a  (3)  was  less  than 
3  percent. 

The  computer  control-program  for  the  spectrometer  was  designed  to  be 
flexible  in  order  to  allow  changes  in  data  collection  procedures.  The  program 
is  written  in  BASIC  with  about  500  lines  of  code  to  control  the  readings  of 
eight  temperatures ,  two  pressures,  and  a/b  values  from  the  waveform  processor 
and  store  the  data  on  a  magnetic  medium  for  future  processing.  The  program  is 
time  controlled.  The  fastest  acquisition  time  for  a  complete  measurement 
cycle  was  30  seconds. 

Measurements  begin  by  starting  the  internal  program  of  the  waveform- 
processor,  which  does  100  averages  on  the  output  A(t)  from  the  InSb 
detector.  Since  averaging  takes  about  14  seconds,  the  computer  program  con¬ 
tinues  data  recovery  from  the  temperature  sensors  which  measure  chamber  (Tc), 
water  bath  (Tw),  room  (Tr),  gas  (Tj  ^  3)*  mirror  (Tm)  temperatures,  and  from 
the  total  pressure  meter  (P).  By  this  time  the  oscilloscope  is  nearly  ready 
with  its  data  so  the  computer  program  returns  to  wait  for  the  final  steps  in 
the  scope  program.  After  the  scope  completes  the  100  averages  on  A(t)  it 
normalizes  the  trace  to  remove  any  baseline  slope,  then  stores  the  waveform. 
Next,  an  auxiliary  output  from  the  scope  switches  the  sweep  generator  voltage 
from  low  (0.36  V)  to  high  (36  V).  The  scope  input  changes  to  the  power  en¬ 
velope  of  the  klystron  with  the  lock-on  disabled.  The  scope  program  performs 
50  averages  on  the  power  curve  and  stores  the  waveform  A^.  The  next  computer 
step  is  to  call  for  the  a/b  and  a^/b^  values  of  the  stored  waveforms. 


All  of  the  measured  data  are  temporarily  stored  in  computer  memory  and 
permanently  stored  along  with  time-of-day  on  magnetic  tape  or  disk  storage  for 
safe  keeping  and  later  recovery  and  processing. 

Real  time  data  were  plotted  on  the  CRT  of  the  computer  to  follow  the 
progress  of  an  experiment.  The  selection  of  what  was  displayed  is  part  of  the 
program  configuration.  Two  typical  examples  of  1-hour  operating  periods  are 
shown  in  Figure  2.  Figure  2a  demonstrates  the  automatic  calibration 
to  or  »  0  dB  for  input  power  variations  of  1  dB.  The  reference  power  trace  aK 
was  reduced  by  a  factor  0.05  for  display,  but  the  correct  coupling  factor 
k  -  0.0550  was  applied  for  the  relative  attenuation  trace  (1) 
a  «  8.686  (/k- a^/a  -  1)  -  0  dB.  Each  1-hour  frame  represents  an  average 
of  8.96  x  10°  actually  acquired  data.  Detection  sensitivity  and  long-term 
stability  are  displayed  in  Figure  2b. 

Program  configuration  is  stored  in  a  separate  file  on  disk.  An  auto¬ 
start  function  is  available  for  studies  that  are  conducted  unattended.  In  the 
event  of  a  power  failure  the  program  will  automatically  re-start. 

Post-processing  of  data  repeats  the  calibrated  signal  correction  as  done 
in  the  real-time  mode.  Scatter  plots  of  ap  are  made  versus  pressure  in  a 
point-mode. 

2.2  Moist  Air  Attenuation  Measurements 

The  objective  of  the  experiments  was  to  perform  pressure  scans  of  the 
attenuation  rate  a  due  to  water  vapor  absorption  in  moist  air.  An  extensive 
series  of  controlled  measurements  was  performed  at  137.8  GHz  to  determine 
a(P  -  e+p).  Over  2000  individual  data  points  a(T,  e,  p)  were  taken  covering 
the  following  range  of  parameters: 

temperature  T  -  316  to  282  K; 

vapor  pressure  e  -  0  to  ej  (RH  >  95%);  and 

total  pressure  P  -  +  p,  where  p  -  0  to  1200  torr  (capacitance 

manometer),  p  -  0  to  3  atm  (aneriod  manometer)*. 
Maximum  vapor  pressure  e^  was  determined  by  the  temperature  Tw  of  the  water 
reservoir  [see  (2 ) -A] . 


♦Experimental  pressure  scale  is  measured  in  torr,  the  prediction  model  MPM 
uses  the  pressure  unit  1  kPa  -  10  mb  -  7.5006  torr. 
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Figure  2.  ftn-wave  spectrometer  operational  data  (5  temperatures  T;  peak 

values  of  reference  and  reflected  signals  a^,  and  a;  and  relative 
attenuation  a,.)  for  a  one-hour  measurement  period  with  an  evacuated 
cell  (fR  -  137.8  GHz):  (a)  automatic  calibration  to  04.  -  0  dB  for 
changes,  Aa*  -  0-1  dB;  (b)  detection  stability  and  sensitivity 
(a  1  -  O.OOS/O.13  <  0.05  dB/km. ) 
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With  the  spectrometer  performance  optimized  at  P  -  0  (see  Figure  2b),  an 
additional  set  of  problems  appeared  when  the  gas  pressure  was  varied.  Intro¬ 
ducing  and  removing  gas  from  an  enclosure  changes  the  temperature  Tg  of  a 
sample  (Figure  3).  Only  pressure  scan  rates  below  ±  100  torr/min  ensured 
quasi-static  gas  conditions,  Tg  -  T^.  Typically,  the  pressure  was  varied  in 
steps.  While  the  gas  settled,  the  klystron  frequency  f^  was  retuned  to 
balance  the  base-line  of  A(t). 

Working  with  water  vapor  often  brought  disappointing  results  with  respect 
to  reproductibil ity .  One  source  was  the  "piston"  effect  where  local  compres¬ 
sion  condenses  part  of  the  vapor;  another  error  source  was  the  slow  diffusion¬ 
mixing  of  water  vapor  with  stagnant  air.  We  calculated  the  diffusion  time 
constant  for  vapor  molecules  to  travel  30  cm  inside  the  cell  against  1  atm  of 
dry  air  to  be: 

T  (  K  )  315  300  285 

Tg  (min)  5.4  6.0  6.7  . 

It  takes  a  period  longer  than  5*tg  for  a  homogeneous  moist  air  mixture 
P  -  ej  +  p  to  develop.  A  measurement  of  a(P)  shown  in  Figure  4a  indicates 
even  longer  time  periods.  Water  pressure  ej(T^)  settled  with  no  delay  when 
the  H20-valve  was  opened.  Dry  air  injection  first  reduced  ej  (piston  effect) 
and  then  it  took  up  to  1  hour  to  obtain  a  stable  value  a(ej+p).  Mixing  was 
accelerated  to  less  than  5  minutes  by  installing  a  fan,  driven  by  a  magneti¬ 
cally-coupled  rotary  feed-through. 

One  other  effect  was  observed  when  dealing  with  moist  air  inside  a 
vessel:  reducing  the  total  pressure  P  caused  the  water  vapor  to  condense  by 
decompression  cooling.  Even  with  the  water  vapor  supply  fully  exposed  to  the 
air  mixture,  it  took  a  long  time  (20  min)  to  reestablish  the  initial  vapor 
pressure  e^  as  indicated  by  ar(P)  in  Figure  4b  (the  mixing  fan  was  off). 
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Figure  3.  Temperature  responses  inside  the  spectrometer  cell: 

(a)  time-response  test  of  gas  temperature  (TQ)  sensors, 

(b)  temperature-vs. -pressure  behavior  for  a  typical  dry  air  case. 
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Figure  4.  Mixing  effects  of  H20  *  AIR  from  time  series  of  138  GHz 

$  decompression  ISldeis^lorduH^iump^iown  over  open  vapor 
source. 


2.3  Water  Vapor  and  Moist  Air  Attenuation  Results 

Moist  air  attenuation  a  at  a  frequency  f  that  falls  within  the  milli¬ 
meter-wave  window  range  centered  at  140  GHz,  can  be  expressed  by  [14],  [17] 

a  -  ks(T)  e2  +  kf (T)  ep  +  kd(T)  p2  dB/km,  (5) 

where  e  and  p  in  kPa  are  partial  pressures  of  water  vapor  and  dry  air,  respec¬ 
tively.  Pressure-broadening  theory  of  the  H2O  rotational  spectrum  predicts 
(e  >  0.01)  a  fixed  ratio  m  between  air-(ep)  and  self-(e2)  broadening;  i.e., 

m  *  kf/kg.  (6) 

An  extensive  series  of  controlled  laboratory  measurements  was  performed  at 
137.8  GHz  to  determine  the  k-coefficients  of  (5)  and  (6).  Table  I  is  a  sum¬ 
mary  of  over  2500  individual  data  points  a(T,  e,  p). 

At  T  -  303  K,  the  foreign-gas  broadener  AIR  (p)  was  replaced  by  its 
principal  constituents  N2,  02,  and  Ar.  These  results  are  summarized  in 
Table  II.  The  broadening  efficiency  m  is  useful  to  explain  H20  absorption 
processes  that  support  (5)  since  a  -  k$e  (e+mp).  Pure  oxygen  (02)  measure¬ 
ments  for  pressures  up  to  2.4  atm  against  Ar  as  the  "loss-free"  reference 
provided  an  estimate  of  kd  when  multiplying  the  result  by  0.21. 

Oata  of  Table  I  were  further  reduced  to  a  reference  temperature 
T0  -  300  K  or  inverse  T-parameter,  0  -  300/T.  The  temperature  dependence  of 
the  ks,f,d  coefficients  was  fitted  to  a  power  law, 

k(T)  -  k0x  dB/km-kPa2.  (7) 

The  results  for  moist  air  attenuation  at  137.8  GHz,  vhen  expressed  by  (5) 
to  (7),  led  to 


x  s:  x  2:  x  x  x  5:  x  3: 


TABLE  I. 


Comparison  Between  Measured  (X)  and  Model- predicted  (M-MPM,  see  Section  3) 
Coefficients  kS(f  of  (5).  Experimental  Conditions:  f*137.8  GHz,  T-282-316  K, 
P-ej  ♦  p,  p-0-1 10  kPa. 


T 

e 

,(RH) 

ks 

Moist 

kf 

Air 

m 

Dry  Air 
kd 

K 

kPa 

dB/km- 

kPa* 

dB/km-kPa* 

X10 

-2 

x10“a 

x10“* 

315.5 

7.49 

(90XRH) 

8.01 

0.485 

6.06 

- 

7.85 

0.481 

6.13 

1.93 

305.9 

4.45 

(90XRH) 

10.9 

0.540 

4.95 

- 

10.81 

0.530 

4.90 

2.10 

303.2 

3.80 

(90XRH) 

12.0 

0.558 

4.65 

2.2 

1 1 .84 

0.545 

4.60 

2.1  1 

296.1 

2.51 

(90XRH) 

15.0 

0.59 

3.9 

- 

15.08 

0.589 

3.91 

2.29 

286.7 

1.39 

(90XRH) 

21  .0 

0.65 

3.1 

- 

21  .22 

0.649 

3.06 

2.46 

28  .8 

1 .05 

(94XRH) 

25.7 

0.68 

2.65 

- 

25.49 

0.687 

2.70 

2.64 

(8) 


k  (T)  -  O.133(4)01O,3(3) , 
kf(T)  -  5.68(5) 1O-303*O(4) , 
kd(T)  -  2.2(5) 1O“602*8, 

and 

m  -  O.O427/07,3. 

Digits  in  parenthesis  give  the  standard  deviation  from  the  mean  in  terms  of 
the  final  listed  digits.  Typical  examples  of  data  plots  ar(e)  and 
ar(ei  +  p)  are  exhibited  in  Figure  5.  All  experimental  results  supported  the 
formulation  in  equation  (5).  Model  predictions  of  the  experimental  data  are 
given  in  Figure  14  (Section  3.4). 
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TABLE  II. 


Attenuation  Measurements  of  Water-vapor/Air-constituent  Mixtures  (f-137.8  GHz, 
T-303  K,  et-3.80  and  p-0-100  kPa)  Expressed  with  kx-Coeff icients  of  (5),  and 
Corresponding  Broadening  Efficiencies  mx  (6).  Also  given  are  Line-core 
Measurements  m  and  their  Resultant  Predictions  m,  and  k,  for  138  GHz. 


LINE  CORE  ( 300K) 
mi. 


I*  ,  tT  T*  V  .  V 


3.  ATMOSPHERIC  PROPAGATION  MODEL  MPM 


(see  Appendix  for  Details) 

Dry  air  and  atmospheric  water  vapor  are  major  millimeter-wave  absorbers; 
so  are  suspended  droplets  (haze,  fog,  cloud)  and  precipitating  water  drops 
that  emanate  from  the  vapor  phase.  A  practical  model  (designated  program 
code:  MPM)  was  formulated  that  simulates  the  refractive  index  n_  -  n'  -  jn"  of 
the  atmospheric  propagation  medium  for  frequencies  up  to  1000  GHz  [3],  [8], 
[9].  Since  the  interaction  with  a  neutral  atmosphere  is  relatively  weak,  the 
refractive  index  is  converted  into  a  refractivity  in  units  of  parts  per 
million, 

•  (n_  -  1)  10®  ppm. 

3.1  Features  of  the  Program 

A  user-friendly  parametric  program  was  developed  that  calculates  the 
values  of  the  complex  refractivity  for  atmospheric  conditions  as  a  function 
of  the  variables  f,  P,  T,  RH,  wA  (A/B/C/D),  w,  and  R,  as  listed  in  tne  scheme 
below: 
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The  output  of  MPM  are  three  radio  path-specific  quantities, 

attenuation  a(f)  dB/km, 

and 

refractive  delay  0Q  ns/km, 

dispersive  delay  B(f)  ps/km. 

The  height  range  30  to  100  km  is  treated  approximately  excluding  the 
detailed  02-Zeeman  effect  and  trace  gas  spectra  (e.g.,  O3,  CO,  N20,  etc.). 
Program  MPM  is  written  with  extensive  comments  to  run  on  IBM-XT/AT  +  8087 
Coprocessor,  or  equivalent,  microcomputers  [9]. 

Extensive  revisions  of  the  MPM  code  have  been  made  over  the  last  year. 

The  newest  version  is  described  in  the  Appendix.  Details  include  (a)  a 
revised  set  of  overlap  coefficients  a^  £(12) -A]  in  the  dry  air  line  table  for 
the  60  GHz  oxygen  band  based  on  a  molecular  fitting  scheme  reported  by  Rozen- 
kranz  [23];  (b)  a  relative  humidity-driven  growth  model  for  hygroscopic  aersol 
addressing  haze  conditions  for  RH  -  90  to  99.9%  [ ( 3 ) -A] ;  (c)  a  new  double- 
Debye  model  describing  the  complex  permittivity  e(f,T)  of  water  for  f  -  0  - 
1000  GHz,  T  -  -10  to  30°C  [(17)-A];  and  (d)  an  improved  fog/cloud  model  that 
includes  a  dispersive  delay  term  [ ( 16 ) -A] . 

3.2  MPM  Cal ibration 

The  experimental  result  (5),  (8)  at  138  GHz  contains  foremost  contribu¬ 
tions  from  water  vapor  continuum  absorption.  Equation  (5)  was  used  to 
"calibrate"  the  program  MPM  by  enforcing  the  agreement  between  experimental 
and  predicted  data.  The  resultant  continuum,  assuming  an  f2  dependence,  is 
defined  by  (see  (14a)  of  [3]) 

ac  ■  0.1820  f  Ng(f) ,  (9) 

where  the  imaginary  part  of  refractivity  is  given  by  (see  Table  III) 

Ng  -  f (3.57  e2010,5  +  0.113  e  p  03)1O“5  (10) 

and  f  is  in  GHz.  The  comparison  in  Table  I  shows  that  at  138  GHz,  within 
experimental  uncertainties,  a  good  fit  was  obtained. 


In  addition  to  (10),  new  width  data  [15]  for  the  183  GHz  H2O  line  (63 
listed  in  Table  1  of  [3]  was  increased  by  11.6%)  *ere  used  to  update  the  MPM 
program;  also,  b3  parameters  of  the  325  and  380  GHz  H2O  lines  were  increased 
by  10  percent.  Other  MPM  modifications  were  a  change  in  the  nonresonant  line 
width  yq  of  dry  air  from  5.6  x  10"^  to  4.8  x  10“^  [23],  the  elimination  of  the 
roll-off  term  1 /[ 1  +  (f/60)^]  (see  equation  13  in  [3]),  and  the  imposing  of  a 
high  frequency  cut-off  F"(f)  -  0  for  f  2  +  40y  [22]. 

3.3  Interpretation  of  H?Q  Continuum  Absorption 

Since  MPM  employs  a  limited  (s  1  THz)  H2O  line  base,  it  is  of  interest  to 
find  out  which  share  of  ks  and  kf  (8)  can  be  attributed  to  far-wing  behavior 
of  the  rotational  H20  spectrum  extending  beyond  1  THz  [7],  [14],  Results  in 
Table  II  indicate  broadening  efficiency  ratios  n  at  138  GHz;  e.g.. 


q  -  m^O  +  ^J/m^O  +  Ar),  (11) 

which  are  similar  to  those  observed  at  cores  of  the  22  and  183  GHz  lines. 

Consequently,  kf  can  be  interpreted  as  a  foreign-gas  broadening  effect  whereby 

MPM  lines  account  for  about  30  percent  of  the  rotational  H20  spectrum. 

Table  III  lists  the  assessment  at  300  K  (0  -  1).  When  the  line  core 

0  5 

argument  (mL  -  0.208/0  *  )  is  extended  to  ks,  then  a  substantial  share, 

<“1  ■  \  '  v 

a.  =  0.103  e2912-5  dB/km  (12) 

12  5 

is  not  supported  by  the  H2O  monomer  line  spectrum  (0  *  -  4.3%/K). 

At  this  point  one  might  speculate  about  ^0)2  dimer  absorption.  An 
estimate  of  the  partial  dimer  vapor  pressure 

eQ  ~  3.12  x  10'4  e205  kPa,  (13) 

was  obtained  by  fitting  data  on  physical  dimer  properties  given  in  Table  IV. 

A  strong  138  GHz  attenuation  rate, 

a.  ~  330  eD07,5  dB/km,  (I4) 

results  when  (12)  and  (13)  are  combined  (0^  ~  2.5%K). 


TABLE  III. 


Summary  Data  of  H20  Line  Spectrum  and  H20  Excess  Attenuation  at  137.9  GHz, 
T-300  K. 


a)  Experiment-versus-MPM  Predictions 

MOIST  AIR 

kg  x3  kf  Xf  B 


DRY  AIR 
kd  xd 


b)  Local  H20  Lines  in  MPM  (31>  of  c): 

0.802  3.5  0.175  3.8  21.8  0.3 

c)  Complete  H20  Line  Spectrum  (m  -  0.022  assumed  -  see  Table  II) 

2.59  3.5  I  0.57  3.0  I  22  -0.5  I 


d)  Excess  H20  Absorption: 


TABLE  IV. 

Reported  Data  on  (H20)2  Dimer  Concentration  e^/e  Over  the  Temperature  Range 
from  300  to  386  K. 
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3.4  MPM  Predictions 

Features  of  the  user-fr i end  1 y  atmospheric  propagation  model  MPM-N  have 
been  discussed  in  Section  3.1  The  desk  computer  version  is  written  in  IBM 
Professional  FORTRAN  with  extensive  comments  that  guide  the  user  to  appro¬ 
priate  references  for  specific  formulations.  Three  parametric  presentations 
have  been  found  useful  in  practical  applications,  which  are  addressed  in 
separate  subprograms: 

Profiles  Variable  Parameters 


s 


A. 

Frequency 

f 

P.  T, 

RH  ( v )  ,  wA,  w,  R 

B. 

Humid ity 

RH,v 

f,  P,  T, 

,  wA,  w 

C. 

Pressure 

P 

f,  ",  T, 

RH  ( e ) 

v! 

V 


i 


P 


? 
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A  copy  of  the  opening  lines  of  program  MPM-N/A  (frequency  profiles)  is  shown 
in  Figure  6.  The  detailed  structure  of  MPM  comes  best  to  light  in  graphical 
examples.  Typical  sea  level  behavior  of  MPM-pred icted  rates,  a  and  8,  is 
illustrated  in  Figures  7  to  13.  Across  the  millimeter-wave  spectrum 
(Figures  7  and  8)  one  recognizes  more  or  less  transparent  window  ranges  (W1  to 
W5/  separated  by  molecular  resonance  peaks.  Calculations  of  total  radio  path 
attenuation  A[dB],  delay  B[ps],  and  atmospheric  noise  temperature  TA[K] 
require  distributions  of  P(x),  T(x)  along  the  propagation  direction  x  to  be 
known  [6]. 

Spectroscopic  properties  of  the  a ir-broadened  water  vapor  line  centered 
at  183  GHz  are  demonstrated  in  Figure  9.  Accuracy  of  supporting  line  para¬ 
meters  (stength,  width,  shift)  determines  the  reliability  of  predictions 
needed  for  remote  sensing  applications  with  respect  to  atmospheric  humidity. 
With  respect  to  line  shape  studies  one  observes  that  wing  data  from  the  delay 
spectrum  8(f)  are  unaffected  by  dry  air  pressure  p.  Outside  a  line-center 
frequency  range  vQ  ±  5y,  dispersion  N * ( f )  is  independent  of  a  particular  shape 
function  F’(f)  such  as  (9b)-A. 
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Modeling  Atmospheric  Millimeter-wave  Properties 

PROGRAM  MPM-N 
JULY  1987 


COMPLEX  RADIO  REFRACTIVITY  OF  ATMOSPHERIC  AIR  (1  TO  1000  QHz). 
Han*  J.  Llebe 

Institute  for  Telecommunication  Sclancaa 

NTIA/ITS.S3 

325  BROADWAY 

Boulder,  CO  80303 

Adapted  for  IBM-PC  Profeaaional  FORTRAN  by: 

John  Strlcklen  (1-303-497-3195,  FTS  320-3195) 

Contents : 

A.  FREQUENCY  PROFILES 
£.  PRESSURES  PROFILES  OF  MOIST  AIR 
S.  HUMIDITY  PROFILES 


A.  FREQUENCY  PROFILES 


Program  MPM-N/A  calculates  fraquency  profiles  of  the  complex  radio 
ref ractivity  N  of  ATMOSPHERIC  AIR  over  a  frequency  range  f = 1  to 
1000  QHz.  The  output  is  expressed  in  Real  part  of  N  and 
Imaginary  part  of  N  and  specific  ratee  of  power  attenuation  A 
(Im.  of  N)  and  propagation  delay  B  (Re.  of  N). 


OUTPUT : 


* 

Real  part  of  refractlvity 

N’  (f) 

in 

ppm 

s 

Imaginary  part  of  rafractivity 

N”(f) 

in 

ppm 

* 

Nondlspersive  Rafractivity 

No 

in 

ppm 

* 

Attenuation 

A(f) 

in 

dB/km 

* 

Dispersive  Delay- 

B(f ) 

in 

pa/km 

* 

Refractive  Delay 

BO 

in 

pa/ km 

Frequency  range,  moist  air,  rain  conditions  (specified  by  five 


meterological  variables  P=p+e 

T,  RH, 

wA  or  w, 

RR),  and  haze  model 

cod#  (A, 

,  B,  C,  or  0)  are  the 

required 

input  information  to  be 

•nt#r#d 

* 

from  the  keyboard. 

INPUT: 

Frequency  range 

f  1 ,  f  2 

in  QHz 

valid  Range: 

1  to  1000 

and  step  size 

df 

in  QHZ 

>=  l.E-05 

* 

Barometric  preasur* 

P 

in  kPa 

(max.  500  freqa. ) 
0.0  to  120 

* 

Temperature 

T 

in  C 

-50  to  50 

* 

Relative  humidity 

RH 

in  X 

0  to  100 

[absolute  humidity 

is  calculated  as 

v ( RH , T )  in  g/m3 

* 

Haze  model  (RH  =  80 

to  99. 9*) 

or 

* ( RH , T  )  in  kPa. ] 

code  A( rural),  B( urban),  C( maritime ) ,  and  D(C  ♦strong  wind) 
plus  hygroscopic  aarosol  reference 
concentration  wA(SOXRH)  in  mg/m3  0  to  i 

«  Suspended  water  droplet 

concentration  w  in  g/m3  0  to  10 

«  Rain  Rat*  RR  in  mm/hr  0  to  200 

COMMENTS: 

The  computation  of  attenuation  A  and  delay  B  is  described  in 
References  [i]  to  [3]. 

C  MOLECULAR  affects  due  to  oxygen,  nitrogen,  and  water  vspor 

Figure  6,  Opening  lines  {total  -  1350)  of  MPM-N/A  PC-computer  program 
listing. 
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Figure  9.  Pressure-broadening  (AIR)  example  of  the  183  GHz  H?0  line. 
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A  major  concern  for  most  NWW  systems  is  their  performance  in  rain.  A 
simplified  classification  of  rain  events  is  given  in  Section  1.1  of  the  Appen¬ 
dix.  When  neglecting  any  statistical  nature  of  rain  within  a  radio  path,  the 
calculation  scheme  ( 18) -A  provides  a  rough  estimate  on  a(f)  and  8(f).  Pre¬ 
dictions  (Figure  10)  are  based  on  adding  to  the  known  state  of  moist  air  (P-T- 
RH)  only  one  additional  parameter,  namely  the  point  rain  fall  rate  R. 

Of  some  importance  is  the  fog/cloud  prediction  program  in  MPM.  Systems 
operating  in  the  NMMW  range  offer  an  attractive  alternative  to  electro-optical 
systems  when  operation  has  to  be  assured  during  periods  of  adverse  weather 
(rain,  cloud,  fog,  haze,  high  humidity).  Current  interest  is  focused  on  the 
frequency  range  90  to  350  GHz  vrfiere  an  optimum  trade-off  between  atmospheric 
obscurations  and  angular  resolution  can  be  achieved  [11].  All  atmospheric 
loss  and  delay  effects  have  to  be  known  accurately  in  order  to  analyze  the 
potential  for  all-weather  performance.  The  suspended  water  droplet  (SWD) 
formulation  ( 16 ) -A,  (17) -A  of  MPM  is  an  addition  to  the  state  of  saturated  air 
( P-T-RH-100%)  as  illustrated  in  Figure  11.  Key  parameters  are  SWD  content  w 
in  g/m3  and  temperature  T  in  °C. 

Another  atmospheric  ingredient  is  the  mass  concentration  wA  of  hygro¬ 
scopic  aerosol  (HAE).  Solution  droplets  appear  for  RH  >  80%,  and  haze  condi¬ 
tions  develop  as  RH  approaches  100  percent.  These  conditions  can  be  modeled 
by  assuming  that  at  the  reference  humidity,  RH  -  80%,  the  HAE  concentration 
wa(RH  -  80%)  is  known.  The  RH  dependent  swelling/shrinking  w(RH)  is  described 
approximately  up  to  RH  -  99.9%  by  a  growth  model  (3 ) -A. 

Relative  humidity  RH  is  the  variable  that  governs  physical  processes 
taking  place  in  the  atmosphere  with  respect  to  water  vapor.  Practical  models 
for  RH  parameterization  exceeding  saturation  are  not  available.  Both, 
absolute  humidity  (v),  and  suspended  haze  droplet  concentration  (wA),  hov«ver, 
can  be  expressed  in  terms  of  RH  variability.  Absolute  (v)  and  relative  (RH) 
humidity  are  interrelated  through  the  ambient  temperature  T. 
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Figure  10. 


Attenuation  a  and  delay  8  rates  for  three  rain  events  (R-10,  50, 
and  100  mm/h)  added  to  a  moist  air  (RH-90%)  sea  level  condition. 
Also  shown  are  dry  air  (OX  RH)  and  moist  air  (50%  RH) 
characteristics. 
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Sulfates  are  the  major  RH-active  ingredient  of  both  urban  and  rural 
aerosols;  so  is  sodium  chloride  for  maritime  species.  The  atmosphere  is  never 
free  from  HAE,  with  greatest  concentrations  near  the  surface  and  scale  heights 
on  the  order  of  1  km.  Above  RH  -  90%,  suspended  water  droplets  have  developed 
carrying  the  HAE  essence  in  solutions.  Average  values  of  wA  lie  between  0.01 
and  1  mg/m^.  The  humidity  parameterization  in  MPM  is  demonstrated  by  the 
examples  given  in  Figures  12  and  13. 

Pressure  variability  comes  into  play  when  modeling  height  dependencies. 
Cumulative  calculations  of  a/0  for  a  slanted  radio  path  through  the  neutral 
atmosphere  (e.g.,  ground-to-satellite)  encounter  pressures,  P  -  100  to  0  kPa, 
which  narrows  the  molecular  absorption  lines  until  they  vanish  altogether. 
Pressure-,  Zeeman-,  and  Ooppler-broadening  [(11) -A]  have  to  be  considered  over 
the  height  range-.O  to  90  km.  Another  need  for  a  formulation  of  pressure 
profiles  arises  from  spectroscopic  studies  applying  pressure-scanning  tech¬ 
niques.  A  simulation  of  laboratory  measurements  discussed  in  Section  2.3  is 
exhibited  in  Figure  14. 
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Figure  11, 


Attenuation  a  and  delay  6  rates  for  three  fog  events  (w  -  0.10, 
0.25,  and  0.5  g/irr)  added  to  a  saturated  air  (RH-100%)  sea  level 
condition.  Also  shown  are  dry  air  (OX  RH)  and  moist  air  (50%  RH) 
characteristics:  (al  f  -  1-1000  GHz, 
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ATTENUATION,  dll /km  DISPERSIVE  DELAY,  pa/km 


320  GHz 

P=101 . 3  kPa 
wA=0.001  ?/h3 

HAZE  MODEL 

4.  EXPE R IMENTAL- VERSUS-MODEL  (MPM)  DATA 


Corroborative  experimental  data  of  sufficient  quality  to  scrutinize  MPM 
predictions  are  scarce.  Reliability,  precision,  and  limited  scope  of  support¬ 
ing  meteorological  data  often  compromise  the  accuracy  of  results  deduced  from 
field  observations.  Generally,  laboratory  experiments  are  more  accurate  by 
simulating  controlled  electromagnetic  and  atmospheric  conditions  crucial  to 
validate  a  specific  model  aspect.  In  this  manner,  contributions  from  the  22 
and  183  GHz  water  vapor  lines  and  from  the  48  to  70,  119  GHz  oxygen  lines  have 
been  evaluated  [3],  [13],  [15],  [16],  [23].  In  many  cases,  theoretical 
refinements  are  motivated  to  improve  the  interpretation  of  empirical 
laboratory  data  by  establishing  a  connection  to  the  physics  of  the  problem. 

For  example,  a  set  of  our  unpublished  dispersion  (N‘)  results  on  dry  air, 
taken  during  1976  between  53.6  and  63.6  GHz,  provided  reference  data  for  an 
elaborate  reformulation  of  interference  coefficients  that  describe  the  60  GHz 
O2  bandshape  [23].  The  new  coefficients  were  adapted  for  MPM,  and  Table  V 
elucidates  the  degree  of  agreement. 

4.1  Laboratory  Measurements 

The  results  on  moist  air  attenuation  at  138  GHz  [i.e.,  (5),  (8)]  provided 
clues  to  a  formulation  of  an  empirical  MPM  water  vapor  continuum.  Data  from 
other  workers  were  evaluated  to  check  if  the  assumptions  of  (10)  held  up  at 
different  frequencies.  Water  vapor  attenuation  at  frequencies  between  330  and 
430  GHz  is  compared  with  MPM  predictions  in  Figure  15.  The  data  were  avail¬ 
able  in  graphic  form  and  had  to  be  digitized  by  us.  Considering  the  many 
difficulties  that  plague  absolute  calibrations,  a  comparison  of  model-vs.- 
experiment  is  encouraging.  No  anomalous  absorption  features  have  been 
uncovered . 

Two  laboratory  studies  of  the  183  GHz  water  vapor  line  have  been  reported 
by  a  French  group  [15],  [16].  We  adapted  the  width  results  [15]  for  the  line 
table  of  MPM.  A  detailed  analysis  of  temperature-dependent  wing  data  [16]  is 
given  in  Table  VI.  Equation  (5)  was  applied  to  reduce  both,  experimental  and 
MPM  data,  for  a  mutual  comparison.  At  +1  GHz  from  the  line  center,  the 
attenuation  rate  for  pure  water-vapor  displays  a  21  percent  discrepancy.  The 
origin  of  such  disagreement  has  to  be  attributed  to  183  GHz  line  broadening 
problems  but  not  to  the  rather  small  contribution  from  the  continuum. 


TABLE  V. 


Measured  (EXP)  and  Predicted  (MPM)  Dispersion  in  Dry  A  r  at  300  K  (MPM  uses 
interference  coefficients  from  Rosenkranz  [23]). 
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Figure  15.  Moist  air  attenuation  rates  a  across  the  atmospheric  window 
ranges  W5  and  W6  (320  to  430  GHz)  at  temperature  8.5,  25.5,  and 
32.9°C:  measured  data  [25],  [26];  solid  lines,  MPM. 
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4.2  Field  Measurements 

Ultimately,  it  is  up  to  field  experiments  to  garner  realistic  evidence  in 
support  of  model  predictions.  When  conducted  with  care  at  different  locations 
under  a  variety  of  natural  conditions,  such  experiments  prove  to  be  quite 
costly. 

Data  from  ITS  field  studies  on  the  propagation  of  millimeter  waves  were 
evaluated  for  water-vapor  attenuation  rates  at  96.1  GHz.  Comparison  between 
clear  air  data  and  MPM  predictions  shows  close  agreement,  including  for  the 
temperature  dependence,  as  demonstrated  in  Figure  16.  The  mission  of  these 
experiments  was  to  establish  a  data  base  for  11.4  to  96.1  GHz  propagation  in  a 
humid  climate  (Huntsville,  AL:  April  -  August  1986).  A  similar  comparison  was 
not  so  favorable  for  a  set  of  data  displayed  in  Figure  17.  The  results  were 
obtained  over  a  1.5  km  1 ine-of-s ight  path  in  the  Soviet  Union  at  sub-freezing 
temperatures  for  frequencies  between  192  and  260  GHz  [24].  Experimental 
uncertainties  are  not  discussed,  the  extreme  environmental  conditions  may  have 
played  a  role.  On  the  other  hand,  335  -  420  GHz  field  data  (included  in 
Figure  15),  that  were  reported  by  another  group  from  the  Soviet  Union  [25], 
agree  remarkably  well  with  MPM  predictions. 
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Figure  16. 


Moist  air  attenuation  o(v)  measured  at  96.1  GHz  over  a  21.4  km 
line-of-sight  path  located  in  Huntsville,  AL  (h  1  0.3  km)  for  6 
temperature  groups  between  2.5  and  37.5°C  [Manabe  et  al.,  private 
communication,  1987]:  x«  5  min  averages  for  4.5  days  (5/4-6, 
8/15-16/1986);  solid  lines,  MPM. 
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•or  attenuation  rates  o(v)  across  the  atmospheric  window 
at  four  temperatures,  5,  -7,  -10,  and  -18°C: 
its  [24];  solid  lines,  MPM. 
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Figure  18.  Terrestrial  path  attenuation  at  94  GHz  and  0.65  pm  (visible  light) 
under  nonprecipitating  conditions  as  a  function  of  absolute  (v) 
and  relative  (RH)  humidity  [8]:  (a)  data  points,  5  min  averages 

taken  every  30  min  during  a  period  of  4  months;  (b)  MPM  simulation 
of  (a)  for  absolute  humidity;  (c)  MPM  simulation  of  (a)  over  the 
range,  RH-95  to  100%,  including  haze  model  C. 
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Attenuation  data  at  94  GHz  have  been  recorded  in  a  coastal  region  of  the 
Netherlands  [8].  Data,  when  presented  in  Figure  18  versus  absolute  humidity  v 
display  significant  random  excess  attenuation  over  MPM  predictions.  The  same 
data  plotted  versus  RH  places  all  excess  attenuation  at  RH  >  98%.  Haze  and 
fog  conditions  probably  were  present  as  evident  from  optical  transmission 
data.  The  highest  excess  of  0.8  dB/km  requires  w  0.16  g/m^  which  is  typical 
for  heavy  fog. 

Figure  19  presents  atmospheric  noise  temperatures  measured  against  zenith 
at  two  different  sites  simultaneously  at  10,  33,  and  90  GHz  [27].  Predictions 
with  the  MPM  radio  path  program  reproduce  the  frequency  correlations  quite 
correctly  which  tends  to  confirm  the  fc  assumption  made  for  H20  continuum 
absorption  (10).  In  addition,  model  data  set  a  limit  for  total  precipitable 
water  vapor  V[mm]  carried  by  the  air  mass.  Noise  exceeding  this  limit  prob¬ 
ably  originates  from  suspended  droplets.  Data  trends  stemming  from 
differences  in  the  f-dependence  of  moist  air  and  SWD  absorption  support  such 
assumption. 
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5.  CONCLUSIONS 


Radio  properties  of  the  atmosphere  are  both  a  barrier  and  a  boon  to 
application  programs  in  the  100  -  1000  GHz  (NMMW)  spectral  range.  The  first 
part  of  this  report  gave  a  somewhat  detailed  description  of  a  laboratory 
experiment  that  had  to  apply  latest  advances  in  digital  electronics  and  cryo¬ 
genic  detection  for  the  sole  purpose  of  coming  up  with  no  more  than  two 
attenuation  coefficients,  k$  and  kf  (8)  at  138  GHz.  At  first  sight,  these  two 
quantities  describe  the  attenuation  rate  of  moist  air  over  a  rather 
inconspicuous  range  from  0.1  to  10  dB/km;  however,  a  more  detailed  analysis 
revealed:  the  k-formulation  provides  conclusive  evidence  on  temperature  and 
pressure  dependences  of  the  water  vapor  continuum.  Attaching  an  f^  dependence 
to  the  138  GHz  results  turned  a  frequency-limited  (<1  THz),  propagation  model 
into  a  useful  tool  capable  of  operating  in  frequency,  humidity,  and  pressure 
domains  of  the  atmosphere  (see  Section  3.1). 

The  report  discussed  details  of  contributions  to  modeling  atmospheric 
NMMW  properties,  foremost  the  attenuation  calibration  of  MPM  employing  new 
laboratory  results  on  moist  air.  Water  vapor  continuum  absorption  and,  above 
RH  -  90%,  droplet  attenuation  are  both  affected  by  relative  humidity  RH.  For 
haze  formation,  the  range  up  to  RH  -  99.9%  was  modeled.  So  far,  in  fog  and 
cloud  situations,  when  RH  exceeds  100  percent,  values  for  w  cannot  be  model - 
generated  from  the  atmospheric  water  vapor  pool. 

The  code  MPM  was  readily  accepted  by  radio  scientists  and  engineers. 

About  90  requests  for  copies  of  MPM  have  been  honored  since  January  1986. 
Quotes  from  received  comment  were  encouraging:  "MPM  proved  of  considerable 
value,"  -  "documentation  is  about  the  best  I  have  ever  seen,"  -  "technical 
flow,  model  development  and  general  visibility  for  application  efforts  is 
superior,"  etc.  The  Microwave  Group  of  the  International  Radiation  Commission 
largely  adapted  MPM  as  atmospheric  transmission  code  (i.e.,  30  times 
referenced  in  "Report  of  Microwave  Group  on  ITRA,"  March  1987). 


Although  MPM  was  found  capable  of  predicting  atmospheric  NMMW  propagation 
limitations,  several  shortcomings  still  exist.  They  are,  for  example,  the 
missing  confirmation  for  a  physical  basis  of  water  vapor  continuum  absorption 
(e.g.,  [17],  [18])  dominating  transmission  in  atmospheric  window  ranges 
centered  at  90,  140,  220,  340  GHz  and  higher;  needed  measurements  of 
spectroscopic  parameters  ( i . e . ,  shape,  strength,  width  and  shift)  for  spectral 
lines  of  the  main  absorbers  (O2  and  H2O)  over  the  full  atmospheric  temperature 
range  (300-200  K);  and  a  lack  for  reliable  subfreezing  transmission  data  to 
clarify  the  problems  that  are  indicated  by  Figure  17.  Further  parametric 
(frequency,  pressure,  humidity,  temperature,  gas  composition  such  as  H2O  + 

AIR)  studies  are  proposed  to  realize  the  benefits  obtainable  from  the  high- 
humidity  performance  of  the  spectrometer  that  has  been  perfected  overcoming  a 
great  deal  of  difficulties.  Pressure-broadening  of  the  183  GHz  line  uniquely 
identifies  monomer  behavior  (see  Table  VI).  A  comparative  study  of  wing  (fQ 
±6  GHz)  and  far-wing  (220  GHz)  responses  to  T  and  v(RH)  variations  would  allow 
an  apportioning  of  local  line  (known)  and  continuum  (unknown)  contributions. 
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THE  REFRACTIVE  INDEX  OF  THE  NEUTRAL  ATMOSPHERE 
FOR  FREQUENCIES  UP  TO  1000  GHZ 


Hans  J.  Liebe 

National  Telecommunications  and  Information  Administration, 

Institute  for  Telecommunication  Sciences,  Boulder,  CO 

1.  INTRODUCTION 

Dry  air  and  atmospheric  water  vapor  are  major  millimeter-wave  absorbers; 
so  are  suspended  droplets  (haze,  fog,  cloud)  and  precipitating  water  drops 
that  emanate  from  the  vapor  phase.  A  practical  model  (designated  program 
code:  MPM)  was  formulated  that  simulates  the  refractive  index  n_  -  n '  -  jn"  of 
the  atmospheric  propagation  medium  for  frequencies  up  to  1000  GHz 
(Liebe,  1985;  Liebe,  1987).  The  main  purpose  of  the  program  is  to  express  the 
electromagnetic  properties  of  the  neutral  atmosphere  in  terms  of  available 
and/or  measurable  quantities. 

The  free  propagation  of  a  plane  electromagnetic  wave  at  frequency  f  and 
an  initial  field  strength  EQ  in  an  isotropic  gas  medium  over  distance  L  is 
described  by  the  complex  transmission  factor 

x  *  E(L)/Eq  -  exp[-(2nf /c) (n "  +  jn')L] 

where  n'(f)  and  n"(f)  are  frequency-dependent  measures  of  delay  and  loss,  and 
c  is  the  speed  of  light.  Since  the  interaction  with  a  neutral  atmosphere  is 
relatively  weak,  the  refractive  index  is  converted  into  a  refractivity 

N.  «  (n_  -  1)10®  ppm. 


in  units  of  parts  per  million. 

Refractivity  N  for  moist  air  can  be  obtained,  in  principle,  by  a  summa¬ 
tion  over  all  absorption  features  in  a  given  volume  element.  In  practice, 
various  degrees  of  approximations  are  employed  to  reduce  labor  and  computer 
time  required,  since  the  number  of  contributing  spectral  lines  by  the  dominant 
absorbers  (water  vapor  and  oxygen)  and  by  various  trace  gases  (e.g.,  O3) 
exceeds  10,000.  The  model  MPM  consists  of  local  H2O  (30x)  and  O2  (48x)  lines 
below  1  THz  and  an  empirical  approximation  to  the  contributions  by  H20  lines 
above  1  THz.  To  complete  the  model,  continuum  spectra  for  dry  air,  suspended 
water  droplets  (haze,  fog,  cloud),  and  rain  are  added.  The  supporting 


spectroscopic  data  base  contains  more  than  450  coefficients. 


1.1  Features  of  the  Program  MPM 

A  user-friendly  parametric  program  was  developed  (see  Section  2)  that 
calculates  the  values  of  the  complex  refractivity 

N  -  N0  +  N* (f )  ♦  jN"(f)  ppm 

or  corresponding  path-specific  quantities  of  particular  interest  to  radio 
engineering;  i.e., 

refractive  delay  8  ns/km, 

dispersive  delay  8(f)  ps/km,  and 

attenuation  a(f)  dB/km; 


both  sets  for  atmospheric 

conditions  as 

a  function  of  the  variables  f,  P 

RH,  wA  (A/B/C/D),  w,  and 

R,  as  listed  in 

the  scheme  below: 

Variable 

Symbol 

Validity  Range 

Medium 

frequency 

f 

i  1000  GHz 

barometric  pressure 

P 

120-0  kPa 

Moist 

temperature 

T 

-  50-  +  50°C 

Air 

relative  humidity 

RH 

0-100% 

haze  model : 

RH-80-99.9% 

code  A,  B,  C,  0 

(or  combinations  thereof) 

Haze 

plus  hygroscopic  aerosol 

o 

reference  concentration 

wa(80*RH) 

0-1  mg /nr 

suspended  water  droplet 

w(100XRH) 

0-10  g/m3 

Fog, 

concentration 

Cloud 

rain  fall  rate 

R 

0-200  rnn/h 

Ra  in 

The  height  range  30 

to  100  km  is  treated  approximately 

excluding  th 

detailed  O^-Zeeman  effect 

and  trace  gas 

spectra  (e.g.,  O3, 

CO,  N2O,  etc. 

that  are  discussed  in  chapter  1.2. 4. 3.  Program  MPM  is  written  with  extensive 
comments  to  run  on  IBM-XT/AT  +  8087  Coprocessor,  or  equivalent,  micro  com¬ 
puters. 
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1.2  Physical  Variables 

The  purpose  of  this  section  is  to  relate  measurable  variables  to  N  model- 
specific  variables,  whereby  the  state  of  dry  air  is  described  by  a  partial  air 
pressure  p  in  kPa,  a  relative  inverse  temperature  variable  is 
0  -  300/(273.15  +  T) , 
and  relative  humidity  is  given  by 

RH  -  (e/es)100  -  (v/vs)  100  -  41.51  (e/05)lO(9,834  9-10)S  100%.  (1) 

Equation  (1)  governs  physical  processes  taking  place  in  the  atmosphere  with 
respet  to  water  vapor.  Absolute  (v)  and  relative  (RH)  humidity  are  inter¬ 
related  through 

v  -  (RH/5.752) 06 10^ 10-9,834  9)  -  7.271  e  0  g/m3  (2) 


I 


8 


& 


v^iere  e  is  the  partial  water  vapor  pressure  in  kPa  as  part  of  the  total  (baro¬ 
metric)  pressure  P  -  p  +  e. 

Water  vapor  variability  at  sea  level  (e.g.,  P  -  101  kPa,  0  -  1.016  or 
22°C)  is  typically: 


Dry 

Normal 

Humid 

Saturated 

V 

1 

10 

17 

20g/m'1 

RH 

5 

50 

85 

Suspended  hydrometeors  are  described  by  the  liquid  water  concentration  w, 
which  relates  approximately  to  optical  (0.55  pm)  visibility  U(km).  A 
schematic  categorization  can  be  made  by 


Stratus 

Convective 

Haze 

Fog 

Cloud 

Cloud 

w  £ 

10'2 

10"1 

1 

5  g/m^ 

U  i 

1.1 

0.27 

0.07 

0.03  km 
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Another  atmospheric  ingredient  is  hygroscopic  aerosol  with  a  mass  concen¬ 
tration  w A  in*g/m^.  Solution  droplets  appear  for  RH  >  80%,  and  haze  condi¬ 
tions  develop  as  RH  approaches  100  percent.  The  growing  haze  droplets  can 
reach  values  (w  <  0.1  g/m^),  sufficient  to  contribute  to  medium  losses.  Haze 
conditions  are  modeled  by  assuming  that  at  the  reference  humidity,  RH  -  80%, 
the  concentration  wA{RH-80%)  is  known.  Any  RH-dependent  swelling/shrinking  of 
w ( RH)  up  to  RH-99.9%  is  modeled  by 

w  -  wa(C1-RH)/C2(100-RH)  -  wA  g(RH),  (3) 


where  g(RH)  is  a  growth  function.  The  following  values  have  been  reported  for 
Cl,  C2 : 


1  - 

Case 

Aerosol  Species 

Cl 

C2 

g ( RH— 99 . 9%)  | 

1 

A 

Rural 

117 

1.87 

1 

91 

B 

Urban 

128 

2.41 

117  ! 

;  C 

Maritime 

183 

5.31 

162  ! 

D 

1  - 

C  +  Strong  Wind 
(>10  km/h) 

197 

5.83 

167  ! 

1 

Average  values  of  wA  for  given  stabilized  climatic  situations  can  be 
found  in  the  literature.  Typical  values  for  wA  lie  between  0.01  and  0.5  mg/m^ 
and  the  increase  to  w  can  be  substantial  as  indicated  by  the  maximum  prefog 
values  g(99.9%). 

Precipitation  originates  as  a  statistical  event  within  clouds  suspended 
in  saturated  air.  Its  vertical  distribution  is  separated  into  two  regions  by 
the  height  of  the  0°C  isotherm,  which  can  vary  between  6  km  and  ground  level 
depending  on  season  and  latitude.  The  lower  part  is  mostly  liquid  drops,  and 
the  upper  region  consists  of  frozen  particles  with  occasional  supercooled 
droplet-loadings  by  strong  updrafts. 

Point  rain  rates  R  have  proven  useful  in  modeling  rain-induced  N 
effects.  Rain  rate  can  be  related  to  percent  time  tp,  a  given  value  occurs 
over  the  period  of  an  "average"  year;  to  the  effective  rain  cell  extent  LR/L ; 
and  to  the  instaneous  suspended  liquid  water  concentration, 

Wp  -  m-  R  g/m^. 

In  terms  of  these  variables,  a  typical  local  rain  may  be  classified  as  follows 
(horizontal  path,  L  -  10  km): 
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Steady 

Hea  vy 
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R 

1 

5 

20 
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tR 
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0.04 

0.04 

The  simple  coefficient  scheme  reveals  some  fundamentals  of  rain.  Changes 
in  the  factor  m  indicate  rain  rate -dependent  characteristics  of  drop  size 
distributions.  Widespread  steady  rain  occurs  more  uniformly  and  favors  small 
drop  sizes  (Si  mm  diameter)  which  stay  in  the  air  longer.  Heavy  showers  are 
more  localized,  favor  larger  drops  and  occur  less  frequently. 


57 


C:  *' 


2.  MODEL  FOR  COMPLEX  REFRACTIVITY 


The  complex  refractivity  in  N  units  (i.e.,  ppm  *  10"^) 

N(f;P/T/RH,wA/w,R)  -  N0  +  N'(f)  +  j N “ ( f )  ppm  (4) 

is  a  macroscopic  measure  of  interactions  between  radiation  and  absorbers.  The 
refractivity  consists  of  a  frequency-independent  term  N0  plus  various  spectra 
of  refractive  dispersion  N'(f)  and  absorption  N"(f). 

In  radio  engineering  it  is  customary  to  express  the  imaginary  part  of  (4) 
as  power  attenuation  rate  a  and  the  real  part  as  propagation  delay 
rate  3  (with  reference  to  vacuum)  as  follows: 

a  -  0.1820  f  N" (f)  dB/km  (5a) 

and 

6  -  3 . 336[Nq  +  N ' ( f) ]  ps/km,  (5b) 

vrfiere  frequency  f  is  in  gigahertz  (GHz)  throughout. 

Radio  refractivity  is  defined  to  be  N  ■  N0  at  f  *  0  and  consists  of  four 
terms;  i.e., 

No  -  N5  *  Ne  +  Nw  +  nR* 

The  individual  contributions  are  described  for  dry  air  by 

N*  -  2.588  p  0,  (6) 

for  water  vapor  by 

N«  -  2.39  e  8  +  41.6  e  02,  (7) 

for  the  SWD  term  N“  by  (16c),  and  for  the  rain  term  by  (18c).  The  results 
(6)  and  (7)  have  been  determined  experimentally  at  microwave  frequencies  where 
dispersive  contributions  N(f)  are  negligible. 

A  calculation  of  the  spectrum  ^(f)  for  frequencies  up  to  1000  GHz 
consists  of  several  additive  parts: 

•Resonance  information  of  na«48  oxygen  lines  and  n^-30  water  vapor  lines. 
•Nonresonant  02  and  pressure-induced  N2  absorption  (Np). 

•Continuum  absorption  from  far-wing  contributions  of  strong  H20  lines  falling 
in  the  frequency  range  1-30  THz  (N^)* 

•Suspended  water  droplet  SWD  term  (f^). 

•Rain-effect  approximation(Np) . 
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Absorption  and  dispersion  spectra  are  obtained  from  line-by-line  calcula¬ 
tions  plus  various  continuum  spectra  _Np>ejW^R  according  to 


Dry  Air  Water  Vapor  SWD  Rain 


N'(f)  -  £  ( SF ' ) i  ♦  +  £  ( SF ' ) i  ♦  N '  +  N;  +  (8b) 

where  S  is  the  line  strength  in  kilohertz,  and  F‘  and  F"  are  real  and 
imaginary  parts  of  a  line  shape  function  in  GHz-*. 

2.1  Local  Line  Absorption  and  Dispersion 

The  Van  Vleck-Weisskopf  function  is  modified  to  describe,  to  first  order 
line  overlap  effects  (Rosenkranz,  1987),  which  leads  to  local  absorption  and 
dispersion  line  profiles  in  the  form 


with  the  abbreviations  A  -  yf/vo, 

X  -  (v0  -  f)2  +  y2»  Y  -  (vQ  ♦  f)2  ♦  y2.  Z  -  (v2  ♦  y2)/vo  . 

Standard  line  shapes  F" (f),  including  the  modified  Van  Vleck-Weisskopf 
function  (9),  predict  in  frequency  regions  of  local  line  dominance  about  the 
same  results  for  N"(f)  as  long  as  F"(f)  exceeds  by  0.1%  the  peak  value  at 
f  -  vQ.  Far-wing  contributions  of  smaller  magnitude  depend  very  much  upon  the 
chosen  shape  function.  For  f  »,  the  wing  response  of  (9a)  becomes  non¬ 
physical  and  is  cut-off;  i.e.,  F"  -  0  when  (vQ  +  f)  >  40y.  So  far,  no  line 
shape  has  been  confirmed  that  predicts  absorption  intensities  over  ranges  10"^ 


to  <  10“®  of  F"(v0).  Far-wing  contributions  from  strong  infrared  water  vapor 
lines,  where  a(v  )  can  exceed  10^  dB/km,  are  accounted  for  summarily  by  empir 
ical  correction  [see  (15a)]. 

The  line  parameters  are  calculated  by  the  expressions  below: 


Symbol 

02  Lines  in  Air 

H20  Lines  in  Air 

Eq. 

S,  kHz 

atP  0‘3exp[a2(l-0)] 

bje  03,5exp[b2(l-0)] 

(10) 

y,  GHz 

a3(p8{0*8"a4^  +  1 . Ie0) 

b3(p0°*6  +  4.8O601*1) 

(11) 

6 

<*C 

a5p  0 

0 

(12) 

Line  center  frequencies  vQ  and  the  spectroscopic  coefficients  a,(*10-7  Hz/Pa) 
to  ag  and  bj(*10-3  Hz/Pa)  to  b3  for  strength  S,  width  y,  and  overlap  correc¬ 
tion  <5  are  listed  in  the  Line  Data  File  of  MPM  (see  part  3:  CD-ROM). 

For  the  02  lines  in  air,  equations  (9)  to  (12)  are  valid  for  altitudes 
h  £  35km  (p>  0.7  kPa),  where  lines  are  pressure-broadened.  Higher  up,  Zeeman 
splitting  and  Doppler-broadening  of  the  Zeeman  components  must  be  taken  into 
account  (Liebe  and  Gimmestad,  1978).  An  estimate  for  h  >  35  km  is  made  by 
geometrically  adding  to  the  pressure  proportional  width  ya  in  (11)  a  second 
term 

yh  -  |y l  ♦  (25H)2J  °*5  GHz  (11a) 

where  H  is  the  scalar  earth  magnetic  field  strength  in  Tesla,  ranging  from 
0.2  to  0.9  x  10  .  The  02  spectrum  vanishes  around  H  -  90  km. 

For  the  H20  lines  in  air  Doppler-broaden ing  has  to  be  considered  at 
altitudes  above  60  km  (p  <  0.07  kPa).  An  approximation  is  made  by  replacing 
the  width  y^  in  (11)  with 

rh  *  (y£  *  Yq)0-5  GHz  (lib) 

where  y^  -  2.14Vg*  10  /8  is  the  squared  Doppler  width.  In  applica¬ 

tions  requiring  the  detailed  mesospheric  line  shape  it  is  necessary  to  apply 
the  more  correct  Voigt  line  shape. 


2.2  Continuum  Spectra  for  Air 

Continuum  spectra  in  (8)  identify  dry  air  and  water  vapor  terms  Np  +  Ng 
and  must  be  added  to  the  selected  group  of  local  (MPM)  O2  and  H2O  resonance 
lines  described  by  (9).  Continuum  absorption  increases  montonically  with 
frequency. 

The  dry  air  continuum 

NJ(f)  -  f  (2a0{y0[l  +  (f/^)2]}"1  +  app  81'5)  p92  (13a) 

and 


N'p(f)  -  aQ{[l  ♦  (f/Y0)2]_1  *1}  P  02  (13b) 

make  a  small  contribution  at  ground  level  pressures  due  to  the  nonresonant  O2 

spectrum  below  10  GHz,  and  a  pressure-induced  N2  spectrum  that  is  effective 

above  100  GHz.  A  width  parameter  for  the  Debye  spectrum  of  0o  is  formulated 

-3  no 

in  accordance  with  (11)  to  be  y0  -  4.8  x  10  (p  +  1 . le) 0  *  (GHz)  (Rosenkranz, 
1987).  The  continuum  coefficients  are  a0  -  3.07  x  10"4  and 
ap  -  1.40(1-1.2  f1-5^"5)^"10. 

Water  vapor  continuum  absorption  has  been  a  major  source  of  uncertainty 
in  predicting  millimeter  wave  attenuation  rates,  especially  in  the  window 
ranges.  Moist  air  attenuation  a  at  a  frequency  that  falls  within  a  window  can 
be  expressed  by 

a  -  ks(T)e2  ♦  kf(T)e  p  ♦  kd(T)p2  .  (14) 

A  series  of  controlled  laboratory  measurements  was  performed  at  137.8  GHz  to 
determine  the  k-coefficients.  Data  a(T,e,p)  were  taken  covering  the  following 
range  of  parameters: 

temperature  T  -  8  to  43°  C 

vapor  pressure  e  -  0  to  e^RH  &  95%)  and 

total  pressure  P  -  e^  +  p,  where  p  -  0  to  150  kPa. 

Experimental  data  were  reduced  to  a  reference  temperature  T0  -  26.85°  C 
(8  -  1).  The  temperature  dependence  of  each  k-coeff icient  (k$  f  d)  was  fitted 

X  *  * 

to  a  power  law  k(T)  -  k6  .  Moist  air  attenuation  at  fx»  137.8  GHz  behaved  as 
fol  lows : 
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k$  -  0. 133(4)01O*3(3) ,  kf  -  5 .68(5) 1O"303 ,0 (4) ,  krf  -  2(1 ) 1O_603-  (14a) 

Values  in  parentheses  give  the  standard  deviation  from  the  mean  in  terms  of 
final  listed  digits.  The  experimental  results  (14),  containing  foremost 
contributions  from  water  vapor  continuum  absorption,  were  used  to  "calibrate" 
the  program  MPM  by  enforcing  agreement  between  experimental  and  predicted 
data;  further,  an  f^  dependence  was  assumed. 

The  water  vapor  continuum  is  derived  from  fitting  experimental  data  (14) 
in  the  case  of 


N“(f )  -  f(bfp  +  bee)e  03  (15a) 

and  based  on  theoretical  data  in  the  case  of 

N'(f)  -  f2b0e  03  ,  (15b) 

where  (Liebe,  1987) 

bf  -  1.13  x  10"6, 

bg  -  3.57  x  10-V-5,  and 

b„  -  6.47  x  10“6  . 
o 

In  summary; 

Equation  (15)  is  needed  to  supplement  local  line  (MPM)  contributions,  the 
coefficient  bf  is  valid  only  for  the  selected  local  line  base  treated  with 

O 

line  shape  (9),  and  the  strong  self-broadening  component  bge  is  nearly  un¬ 
affected  by  (9).  The  coefficient  bQ  and  both  exponents  in  (15b)  were  obtained 
by  fitting  dispersion  results  of  line-by-line  calculations  for  the  rotational 
H20  spectrum  above  1  THz. 

2.3  Suspended  Water  Droplet  Continuum  (Haze,  Fog,  Cloud) 

Suspended  water  droplets  (SWO)  in  haze,  fog,  or  clouds  are  millimeter 
wave  absorbers.  Their  size  range  of  radii  is  below  50  pm,  which  allows  the 
Rayleigh  approximat ion  of  Mie  scattering  theory  to  be  used  for  calculating 
refractivity  contributions  to  (8)  in  the  form 

NJJ ( f )  -  (9/2)w/e"(l  *  n2),  (16a) 

N;(f)  -  (9/2)w  [l/(e0  ♦  2)  -  n/e"(l  ♦  n2)],  (16b) 
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and 


N"  =  (3/2) w  [1  -  3/(e0  *  2)], 


(16c) 


where  n  -  (2  +  £')/e";  e',  e"  are  real  and  imaginary,  and  eQ  static  parts  of 


the  permittivity  for  water.  The  contribution  of  (16c)  is  added  to  equation 

(6). 

Values  for  the  dielectric  spectra  e(f)  of  water  are  calculated  with  a  new 
double-Debye  model  (Manabe  et  a  1 . ,  1987): 


e'(f)  -  e2  ♦  (eQ  -  Ej)/[1  ♦  (f/fQ)2]  +  (£j  -  e2)/[l  ♦  (f/f$)2], 


(17a) 


e"(f)  -  f{(eQ  -  Cl)/fD[l  +  (f/fQ)2]  ♦  (ej  -  £ 2 ) / f s [ 1  +  (f/fs)2]J,  (17b) 


eQ  -  77.66  +  103.3(0  -  1), 


(17c) 


where  e1  -  5.48,  e2  "  3.51, 


20.09  -  142(0  -  1)  +  294(0  -  lp  GHz,  and 


fs  -  590  -  1500(0  -  1) 


GHz. 


Equation  (17)  is  valid  for  frequencies  up  to  1000  GHz  over  a  temperature  range 
from  -10  to  +30°  C. 


2.4  Rain  Effects 


The  refractivity  of  rain  is  identified  in  (8)  by  _Np  -  Np  +  jNp.  Drop 
diameters  (0.1  -  6  mm)  and  millimeter  wavelengths  are  comparable,  thus  causing 
appreciable  interactions  due  to  Mie  absorption  and  scattering.  Bypassing 
elaborate,  lengthy  Mie  calculations  which  require  drop  shape  and  size  distri¬ 
butions  as  well  as  the  complex  dielectric  properties  of  water  (17),  rain 
refractivity  spectra  are  approximated  via 


where  fR  -  (52  -  0.22R)  GHz.  Frequency-dependent  coefficient  a  and  exponent 
b  were  calculated  using  drops ize  spectra  of  Laws  and  Parsons  and  a  temperature 
of  T  -  0°  C.  A  regression  fit  to  individual  (a.b)-pairs  over  the  frequency 
range  from  1  to  1000  GHz  resulted  in  the  following  calculation  scheme: 
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3.  CONCLUSIONS 


The  parametric  model  MPM  for  atmospheric  refractivity 

N(f,P/T/RH,  wa/w,R) 

was  developed  for  applications  in  areas  such  as  telecommunications,  remote 
sensing,  and  radio  astronomy.  Details  of  its  structure  and  operation  are 
explained  in  the  extensive  COMMENTS  part  of  the  code,  that  is  available  under 
MPM  in  part  3:  CD-ROM.  The  memory  capacity  required  for  MPM  is  355  kbytes. 

The  format  of  the  numerical  print-out  is  demonstrated  by  the  identical 
example  given  in  TABLE  1  for  the  _N  version  [see  (4)]  and  in  TABLE  2  for 
the  ct/8  version  [see  (5)].  A  plotting  system  at  the  user's  choice  (e.g., 

HALO)  can  be  added  to  include  features  such  as  auto  or  manual  scaling, 
multiple  cases  (e.g.,  9  curves  with  up  to  500  points  each),  special  labels, 
etc.  An  example  of  a  graphical  presentation  for  a  sea  level  condition  of 
moist  air  (wA  -  w  -  R  -  0)  exhibits  spectra  at  various  relative  humidities 
(RH  -  0-100%)  for  the  _N  version  in  FIGURE  1  and  for  the  identical  case  as  ot/B 
version  in  FIGURE  2. 
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TABLE  1. 

FREQUENCY  PROFILES  OF  ATMOSPHERIC  COMPLEX  REFRACTIVITY 


INPUT  Valid  Parameter  ranges  indicated  by  [  ] ) : 

RAIN 
RATE , R 
(mm/hr) 
[0-200] 

10.0 

Minimum  Frequency  FI  0.000  (GHz) 

Maximum  Frequency  F2  [ 1 000 . ] 1 000 . 000  (GHz) 

Frequency  Step  [max  500]  dF  100.000  (GHz) 


REL. 

HAZE 

SUSP. 

CASE 

PRES. ,P 

TEMP. , T 

HUM .  ,  RH 

MODEL 

DROP. ,w 

(kPa) 

(C) 

(*) 

(mg/m3) 

( g/m3 ) 

[1-9] 

[0  -110] 

[+/-50] 

[0-100] 

[0-1] 

[0-10] 

1 

101.3 

15.0 

100.0 

:  0.00 

1  .000 

OUTPUT : 

Case  Number:  1  (No  =  351.18  ppm) 


MOIST  AIR  (v=  12.81  g/m3) 

DRY  WATER  HAZE,  FOG 

AIR  ♦  VAPOR  +  CLOUD  +  RAIN  =  TOTAL 


FREQUENCY 

(GHZ) 

0.000 

1 00 . 000 

200.000 

300.000 

400.000 

500.000 

600.000 

700.000 

800.000 

900.000 

1000.000 


N'-IMAGINARY  PART  (ppm) 

N'-REAL  PART  (ppm) 

+♦♦+♦+++♦  ♦++♦♦♦♦+♦  ♦++++++++  +♦+♦♦+++♦  +++++++++ 

0 . 220E-09  - . 31 7E-07  0.278E-08  Q.431E-08  -.244E-07 

0 . 1 68E-02  0 . 454E-01  0.242E+00  0.317E+00  0.607E+00 

-.21 9E+00  0.321E+00  -.139E+00  -.226E+00  -.264E+00 

0 . 476E-03  0.149E+00  0.288E+00  0.190E+00  0.629E+00 

- . 1 70E+00  0 . 1 07E+01  -.299E+00  -.266E+00  0.334E+00 

0.561E-03  0 . 1 74E+00  0.284E+00  0.120E+00  0.579E+00 

- . 1 62E+00  0.352E+01  -.391E+00  -.275E+00  0.269E+01 

0 . 807E-03  0.481E+00  0.276E+00  0.889E-01  0.845E+00 

-.157E+00  0.61 9E+01  -.450E+00  -.278E+00  0.530E+01 

0 . 1 04E-02  0.11 8E+01  0.268E+00  0.675E-01  0.152E+01 

-  .  1 62E+00  0 . 1  92E+02  -.495E+00  -.279E+00  0.183E+02 

0.846E-03  0 . 225E+01  0.260E+00  0.549E-01  0.257E+01 

- . 1 59E+00  -.172E+02  -.532E+00  -.280E+00  -.182E+02 

0 . 977E-03  0.101E+01  0.252E+00  0.461E-01  0.131E+01 

- . 1 57E+00  0 . 687E+01  -.563E+00  -.281E+00  0.587E+01 

0.117E-02  0.11 2E+01  0.242E+00  0.396E-01  0.140E+01 

-.161E+00  - . 635E+0 1  -.589E+00  -.281E+00  -.738E+01 

0.101E-02  0 . 795E+00  0.232E+00  0.347E-01  0.106E+01 

- . 1 59E+00  0.716E+01  -.611E+00  -.281E+00  0.611E+01 

0 . 1 03E-02  0 . 603E+0 1  0.223E+00  0.308E-01  0.628E+01 

- . 1 59E+00  - . 1 34E+02  -.629E+00  - . 282E+00  -.145E+02 
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TABLE  2. 

FREQUENCY  PROFILES  OF  ATTENUATION  AND  DELAY  RATES 


INPUT  Valid  Parameter  ranges  indicated  by  [  ] ) : 


REL. 

HAZE 

SUSP. 

RAIN 

CASE 

PRES. ,P 

TEMP. , T 

HUM. ,RH 

MODEL 

DROP. ,w 

RATE,  R 

(kPa) 

(C) 

(*) 

(mg/m3) 

(9/m3) 

(mm/hr) 

[1-9] 

[0  -110] 

[+/-50] 

[0-100] 

[0-1] 

[0-10] 

[0-200] 

1 

101.3 

15.0 

100.0 

:  0.00 

1  .000 

10.0 

Minimum  Frequency  FI  0.000  (GHz) 
Maximum  Frequency  F2  [1000. 31000.000  (GHz) 
Frequency  Step  [max  500]  dF  100.000  (GHz) 


OUTPUT : 

Case  Number:  1  (Refractive  delay  =  1171.5  ps/km) 


MOIST  AIR  (v=  12.81  g/m3) 

DRY  WATER  HAZE,  FOG 

AIR  +  VAPOR  +  CLOUD  +  RAIN  =  TOTAL 


FREQUENCY  <X  -  ATTENUATION  (dB/km) 


(GHz) 

0.000 

0.00 

0.00 

0.00 

0.00 

100.000 

0.03 

0.83 

-0.73 

1  .07 

200.000 

0.02 

5.44 

-0.57 

3.56 

300.000 

0.03 

9.50 

-0.54 

11.75 

400.000 

0.06 

35.02 

-0.52 

20.64 

500.000 

0.09 

107.25 

-0.54 

64.18 

600.000 

0.09 

246.00 

-0.53 

-57.35 

700.000 

0.12 

128.84 

-0.52 

22.92 

800.000 

0.17 

162.65 

-0.54 

-21 .20 

900.000 

0.17 

130.27 

-0.53 

23.88 

1000.000 

0. 19 

1097.36 

-0.53 

-44.84 

DISPERSIVE  DELAY  (ps/km) 


0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

4.41 

5.78 

11.05 

-0.46 

-0.75 

-0.88 

10.50 

6.93 

22.88 

-1  .00 

-0.89 

1.11 

15.52 

6.57 

31 .82 

-1.30 

-0.92 

8.99 

20.11 

6.32 

61 .51 

-1.50 

-0.93 

17.68 

24.43 

6.14 

137.92 

-1  .65 

-0.93 

61.06 

28.44 

6.00 

280.53 

-1.77 

-0.94 

-60.59 

32.07 

5.87 

166.91 

-1 .88 

-0.94 

19.58 

35.28 

5.77 

203.87 

-1 .96 

-0.94 

-24.63 

38.08 

5.68 

174.20 

-2.04 

-0.94 

20.37 

40.50 

5.61  1143.65 

-2.10 

-0.94 

-48.40 
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ps/km  ATTENUATION.  dB/km 


Figure  A2.  Moist  air  attenuation  (a)  and  delay  (s)  rates  for  a  sea  level 
condition  (P , T )  and  various  relative  humidities  ( RH )  over  the 
frequency  range  from  1  to  1000  GHz. 
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